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(Anti-)instanton behaviour in Euclidean non-abelian field of the point-like source is studied 
by analyzing the possible (anti-)instanton deformations as resulted from the variations of its 
characteristic parameters. The variational principle for searching such crumpled topological 
t^j- | configurations is formulated and the problem is resolved then algebraically by the Ritz method 

(multipole expansion of the deformation fields). The region of (anti-)instanton parameters rele- 
■ vant to the instanton liquid model is investigated in detail. The approximate efficacious method 

to include such configurations contributing to the appropriate partition function is proposed and 
verified. Evaluating average energy of Euclidean non-abelian point-like source in the instanton 
environment (liquid) is performed within the superposition ansatz for the corresponding solu- 
' tions of the Yang-Mills equations. For the colour singlet dipole this energy escalates linearly 

with the separation increasing and its tension coefficient develops the magnitude commensurate 
with that as inferred from lattice QCD. 
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r-| Considerable interest in classical theory of non-abelian gauge fields has been provoked by the dis- 
covery of instantons when it becomes clear that conceptually low energy QCD shares many common 
Ph| features with condensed matter physics. These pseudo-particle (PP) solutions of Euclidean classical 
J|h I equations of the Yang-Mills theory, having been properly included in the Euclidean functional in- 
tegral, turned out very suggestive in illuminating the nature of QCD vacuum and non-perturbative 
dynamics pQ (the CDG scenario). In particular, it has been realized the QCD ground state is a 
sophisticated strongly interacting system filled with the condensates of quark-antiquark pairs and 
^ \ gluons. Moreover, experimentally observed particles respect only a part of the full symmetry of QCD. 
Then similarly to the condensed matter physics we have chiefly two theoretical approaches to explore 
these systems. One is less theoretical and based on the numerical simulations. Another relies on 
the construction of effective theories for the low energy degrees of freedom governed by studying the 
patterns of breaking of underlying symmetries. The most advanced example is given by the chiral 
symmetry of QCD in the limit of massless quarks. This symmetry is spontaneously broken by the 
condensation of quarks and anti-quarks in the QCD vacuum. Such a phenomenon (together with 
many consequences) is explained to a large extent if one surmises the classical solutions of Yang- 
Mills equations with nontrivial topology ((anti-)instantons with perturbative oscillations about them 
added) are playing dominant role in formation of the QCD vacuum. According to this scenario in the 
low energy region the coupling constant becomes frozen somewhere at the scale of average instanton 
size and the liquid (or gas) of (anti-)instantons randomly oriented in colour space and homogeneously 
distributed over 4-dimensional Euclidean space is an adequate ensemble to saturate the correspond- 
ing functional integral. The quarks are living in this background and its influence after all defines 
the observable quantities such as the gluon and chiral condensates, the dynamically generated quark 
mass the pion decay constant etc. 
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The pressing need in the profound study of quasi-classical approach at calculating the generating 
functional of QCD comes also from the recent developments of investigating a colour superconductiv- 
ity which appears as a new phase of hadronic matter at very high density of quark/baryonic matter 
and low temperatures [2] . The analysis of quark condensate behaviour in the instanton liquid model 
(IL) jH] shows the corresponding phase transition takes place at a surprisingly small value of the 
critical density of quark/baryonic matter (order (or even less) of the quark density in normal nuclei). 
Applying the IL scheme devised in jl] which takes into account the medium collective excitations 
essentially allows to increase the critical density. Moreover, this density of phase transition may 
become so high that the Coulomb fields generated by quarks turn out comparable to the instanton 
ones due to very small interquark distances. Obviously, it puts in doubt the predictions of IL model 
which ignores the influence of external fields on the pseudo-particles and makes IL hardly applicable 
under such conditions. 

Here we focus on studying the (anti-) instanton behaviour in the field of point-like Euclidean colour 
charge. Actually, this problem is one of the clue elements in constructing a more realistic version 
of IL, in particular, at high quark/baryonic densities. Besides, it can shed light on the challenging 
problem of colour field penetration into IL. We argue and deal with the field theory approach in quasi- 
classical approximation aiming to estimate the leading contributions. Apart from scrutinizing the 
terms of interaction between source and pseudo-particle which have been explored in the pioneering 
papers [l] and where its dipole nature has been argued, we are trying to clarify what happens to a 
pseudo-particle itself when an external field is available. Apparently, the problem of investigating 
the deformation fields is of great importance in this context. In spite of the rather impressive record 
of studying the pseudoparticle interactions (see, for example, 0,0) nowadays we are still far from 
declaring this problem resolved and transformed to the practical instrument. We formulate a further 
new approach to study the interactions of point-like solutions of field equations and suggest a new 
efficient way of approximate calculation of the functional integral for interacting pseudo-particles. 
In this way we construct the perturbative approach intrinsically related to the nontrivial topological 
solution itself analyzing the variations of solution parameters. In the case of pseudo-particles these 
configurations retain the bulk topological features and are sensitive to the presence of perturbative 
factors. At first sight, this rather routine excercise exhibits unexpectedly a wealth of possibilities 
and turns out very intersting even out of the IL context. We find out the pseudo-particle could be 
considered as a 'supplier' of various fields, in particular, a scalar field and colour vector and tensor 
fields. We expect the scalar field could be a good pretender to be an interaction carrier in the 
soft momentum region [7] whereas the vector field carrying the colour indecies could be responsible 
for the screening in the instanton medium. The paper is organized as follows. The motivation to 
consider the problem and supposed method of its treatment are discussed in Sec.l The deliberation 
of deformation contributions to the total action of point-like source and instanton is given together 
with an approximate approach to calculate the pseudoparticle deformations in Sec. 2 and the practical 
realization of this suggestion is demonstrated in Sec. 3. In the fourth Section we present the analysis 
of numerical calculations and explore the generalized formulation of searching the deformations with 
imposed constraints. And eventually Section 5 is devoted to an analysis of the results obtained and 
their applications for calculating within the IL approach. 

1 Optimizing the instanton configurations 

Let us suppose the presence of immovable Euclidean colour point-like source of external field at the 
point z e . Its intensity is supposed to be e and an orientation along the third axis of colour (isotopic) 
space because for the sake of simplicity we limit ourselves here to dealing with SU(2) group only. 
Then the current density takes the form j a = (0, e 5 a3 5(x — z e )) where the 3-vectors 

H = 1,2,3,4 (as other Greek indecies) and a = 1,2,3 (as other Latin indecies). Fixing the location 
of point-like sources is gauge invariant unlike specifing their orientation in isotopic space. In order to 
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avoid such a stain jH] the static solutions of the Yang-Mills equations should be characterized by their 
energy and 'total' isospin because these quantities are gauge invariant in addition to being conserved. 
As known, the field originated by the particle source in 4d-space develops the cone-like shape with 
the edge being situated just at the particle creation. Then getting away from that point the field 
becomes fully developed in the area neighbouring it. In distant area where the field penetrates into 
the vacuum it should be described by the retarded solutions obeying, in particular, the Lorentz gauge. 
Here we are interested in studying the interaction in the neighbouring field area and the cylinder- 
symmetric Coulomb field might be its relevant image in 4d Euclidean space. Certainly, we skip any 
impact of the instanton field on the source of colour field although the self-consistent solution of the 
problem with non-abelian fields teaches about the possible change of the source orientation in colour 
space (see an example of two point-like colour charges iQj), see also ^U]- However, bearing in mind 
that the instanton field is rather short range (the order of the PP size) we believe the source field is 
not reconstructed essentially on that scale (albeit at the end of this paper we will be well prepared 
to formulate the corresponding generalized equations for the field of point-like source). 
The action for gauge field A ^ when an external source is available reads 

S= [ dx (\ G a G a + jtA'l) , (1) 



with the field strength defined in standard way 

G% = MS " d » A l + 9 z abc K A l > ( 2 ) 

where e abc is an entirely antisymmetric tensor and g a coupling constant for non-abelian field. Actu- 
ally, the solution of the Yang-Mills equations corresponding to a single (anti-) instanton, for example, 
in the singular gauge is well defined and characterized by the coordinate of the PP centre z, the PP 
orientation in the colour space uo ab and its size p, i.e. 

A » = \ | , (3) 

where y = x — z and fjj,^ is the 't Hooft symbol (for anti- instanton fj — > rj). When the distance 
between the point-like source and the instanton centre is large comparing to its size A 3> p (A = |A|, 
A = z — z e ) the standard superposition 

a; = b;(x) + a;(x) , (4) 

of point-like source field which for the sake of simplicity is considered in this paper in the simplest 
gauge only 

B;(x) = (0,6°* V ), ^A—I—, (5) 

and the (anti-)instanton field Eq.fjSJ) is a pretty reasonable approximation of the true solution. The 
Lorentz gauge condition d^B^ = is valid for the retarded solution (as for the instanton field). 
In the region of a well developed field this condition is obeyed by the Coulomb solution Eq.© as 
well. We treat the potentials in their Euclidean forms but the following changes of the field and the 
charge of point-like source Bq — > iB^, e —ie should be done at transition from the Minkowski 
space. Actually, last variable change is resulted from the corresponding transformations of spinor 
fields tj) — > if), "tj) — > 70 — > 74 where the hatted spinors are Euclidean. Clearly, these Euclidean 

sources generate the fields of the same nature as ones we face taking into account gluon field quantum 
fluctuations a qu around the classical instanton solutions A = Ai nst + a qu . 

Searching the optimal PP configurations we fix its size and orientation in the isotopic space at a 
large distance from the source. But approaching the source (in the local vicinity of it) both are able 
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to vary as the functions of the PP centre coordinate, i.e. p — *> R(x,z), uo ab — > Q ab (x,z). Now with 
such substitutions in the superposition ansatz Eq.(0} we may formulate the variation problem of 
searching the deformation fields R(x, z), Q(x, z) which optimize the action Eq.(£Q) with the boundary 
conditions for these fields falling down at A — ► oo. Finally, as an output we come to the sophisticated 
system of differential equations which can be analyzed numerically. However, such a scheme even 
treated perturbatively faces the principle difficulty because the coefficients of higher order derivatives 
of the deformation fields R(x, z) and Q(x, z) are strongly suppressed beyond PP owing to the almost 
singular behaviour of the (anti-)instanton solution. This fact essentially complicates the process 
of searching the appropriate solutions and, moreover, makes it practically impossible sometimes. 
At the same time this peculiarity of the PP solution behaviour prompts very natural approach to 
circumvent the difficulty mentioned. It implies applying the variational Ritz procedure (instead of 
the grid method) which in this particular case can be realized as the multipole expansion of the 
(anti-)instanton size, for example, 

R in {x, z) = p + c^y^ + c^ y tl y v + ... , \y\ < L 

Rout{x, z) = p + oL -§ + dn U — 77 — 2 + • • • , \y\ > L . 

V V V 

Similar expressions could be written down for the (anti-)instanton orientation in the colour space 
Q(x, z) where L is a certain parameter fixing the radius of sphere where the multipole expansion 
growing with the distance increasing should be changed for the decreasing one being a result of 
imposed deformation regularity constraint. Hoping to rely further on the perturbative analysis we 
have to calculate the coefficients c M , c^, . . . and cf M , d^, . . . and the parameter L by minimizing the 
corresponding quadratic form resulting from Eq.(JTJ), i.e. 5S = 0, (being calculated up to the second 
order terms) together with the boundary conditions at infinity and the conditions of sewing 'in' and 
'out' parameters on the sphere \y\ = L. The final results allow us to conclude that the instanton field 
turns out to be so strong regularizator that we are allowed to deal with the fields Ri n and Qi n only 
at exploring the terms of the second order in small deviations from p and uj. One may neglect all the 
specifications coming from R out and VL out which means formally to operate in the limit L^> p. If one 
is interested in stationary cylindrically-symmetric picture in the 4-dimensional space the dependence 
on temporal coordinate 24 should be excluded, i.e. the spatial indecies 1, 2, 3 have to be retained in 
the coefficients of multipole expansion Eq.(jH}. Then it is clear the final result becomes dependent on 
the distance between three dimensinal components of the instanton centre and the source position 
which is denoted as A hoping it does not produce a mess. We would like to notice here that our 
condition to defined the deformations 5S = allows us to fix the major contributions to the exponent 
of the generating functional. In a sence we might consider this equation as the valley equation in the 
quasiclassical approximation. Then the higer order terms should reproduce loop corrections (would 
be quantum interactions). 

Going to implement the program mentioned above let us consider the example of the dipole 
expansion terms © denoting them as 

5p = 5^p |/ M , 5u = , R = p + 5p , Q = u + 5uj . (7) 

Then for the crumpled instanton we have instead of (jHJ) the following expansion in 5p up to the 
second order terms 

dAliVt.p) d 2 Al{yi,p) (AR) 2 

A a .. ~ A*i$l, p) + — ^ v HJ AR + , 
n /A >H)-r dR OR OR 2 ' 

here AR = R — p. Using the superposition B a + A a in Eq.(J2J leads to the definition of G a as 



G% = G% / (B) + G%(A) + G%(A 1 B) 
G%{A,B) = g e ahc {B b ^Al + AlBl) , 
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with G a v (A) and G a ^(B) in the standard form of Eq.(j2J) and the field strength squared looks then 
like the following sum 

G% G% = G%,{B) G%(B) + G%{A) G%{A) + G%{A, B) G%(A, B) + 

+ 2 G%(B) G%{A) + 2 G%(B) G%(A, B) + 2 G%(A) G%(A, B) . 9 

One should not envisage this complicated expression tractable in the approximation of weak external 
field because of the strongly singular nature of instanton solution and point-like source field. Besides, 
the large number of interacting terms essentially contributing does not allow one to be based on 
their physical meaning only (as in [TT]) in further analysis. Now let us note the non-zero strength 
component of point-like source field reads 

G ^ B ) = — 5 a3 Vl + Al 



Aix \y + A\ 6 

whereas the field strength generated by the instanton is presented as 

g^a) = g%(A) + ag« u (A) , (10) 

and the first term here looks like instanton field strength but with the varying parameters 

G%(A) = -- n ak M aa M uP f} ka0 (y2 ^ , 

where = 5 ua — 2 , and the correction term is given by 

= dA a v (Q,R) dR _ dA«(Q,R) OR dA a u (Q,R) dn bc _ dA°(Q, R) 3Q bc 
^ [ ' dR dy a dR dy v dtt bc dy„ dtt bc dy v 

+-^e abc (uj hm 5uj cn - uj bn 5u cm + 5cu bm tu cn ) - 5u ak e kmn ) fj m „ a fjn»p f" R 2y • 
Finally, the interfering term takes the following form 

a3c nek = Vol R ' 



G a u {A,B) = 2 — e aAc n ck fj ki 



4tt " wa y 2 (y 2 + R 2 ) \y + A| " 

The other terms do not contribute to G a (A, B) because of the particular choice of gauge. Later we 
will develop the perturbative approach for this part of strength expanding it in 5p 5u up to the 
second order terms. 



2 Classifying the deformation contributions 

The first term of Eq.fjHJ) is certainly out of our interest here and calculating the action with the 
second term included gives the correction to single instanton action. The quadratic form of kinetic 
energy type hinges upon the deformation fields 

i kl 

Skin — j dx — G%(A)G%(A) - /? = | (M 2 + ^ ^ a %u al + ^ e kmn e abc u bm 5,u ak 5,uj cn + 

a a ( n ) 

+y v muu e abc (uj bm 5 u u; cn - co bn 5^ cm )^ an + -j u am u cn 5^ am 6^ cn + v% 5,tu al 5„p, 
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with the tensor coefficients 
9 

K = To ^ ' a ^ = ai 5 ^ 5kl + ° 2 £kln ^ u ' v °» v = ' v % = Vl 5 ^ Skl + V2 £kln ^ ' 

23 3 

^i = ^/V, a 2 = --f3p\ a 3 = -2v lP , (12) 

4 7 1 

cr 4 = 2v 2 p , <x 5 = - v 2 p , «i = — p p , v 2 = - (3 p , 

where (3 = is the single (anti-)instanton action. These contributions are regulated by instanton 

itself and independent of the field B as well as the distance between the pseudo-particle center and 
the source. In a sense they characterize the PP compliantness as for its shape changes and an 
orientation in colour space. Here is an appropriate place to notice that dealing with the crumpled 
configurations characterized by the parameters (J2J) within the perturbative approximation and in 
the limit L — > oo (see (jOJ)) presupposes a lot of singular terms available. Fortunately, the natural 
regularization pops in when those bad terms are multiplied by the instanton tensor G" (p, uj) (see 
below ([22))) . It is interesting to notice the operator similar to Eq.Qlip emerges when one explores the 
fluctuation fields in the proximity of the particular classical solution A = A c i + a qu and the action 
acquires an additional quadratic term 

S = S{A d ) + dx a qu C(A d ) a qu = S(A d ) + ^ A* ^ , (13) 
J k 

where are the eigenvalues of the operator C and presents the coefficients of the expansion of 
field a qu = Y^kikipk into the entire set of orthogonal eigenfunctions ipk of this operator. In view 
of that generally Eq. fjllJ) delineates the parametric dependence of eigenvalues and eigenfunctions of 
the operator £ in the accepted approximation. In our particular case it means dependence on the 

expansion coefficients in Eq.© (|K !*.■■■)■ Apparently, we gain the possibility to take into 

account the contribution of regular component originated by external field in the functional integral 
instead of dealing with poorly defined Green function of instanton. 

Calculating Eq.([T2"|) and some subsequent ones we used the following identities valid for SU(2)- 
group 

uj ab uj ac = S bc , e abc u aa u bp cu^ = e al3 ^ . (14) 

It is practical for analysing the additional term to the single instanton action to present it as a 
sum of two components 

S — (3 = Skin + Si n t ■ (15) 

where Skin is the term of kinetic energy type and Si n t is the interacting term (of course, the singular 
contribution of the point-like source to the action is ignored). The interacting term Si n t may be 
decomposed into the following sum 

Sint = S^t + Sfnt + ^int ■ (16) 

and the physical meaning of each term here becomes clear from further calculations. We quote here 
the final result only ommitting the routine calculations which are not very inventive and require the 
good patience mainly. 

2.1 Quadratic terms depending on the distance from source 

As the first term of Eq. (jl6|) we mean the quadratic terms in Sp, 5u wholesale which describe the 
dependence of kinetic coefficients on A: 

1 e / e \ 2 

Sint = — T~ ^VkH ^lij $iP 5jP + I — (Ski - UJ 3k U 3l )(S k l Qij - fjk4mfjl4n Rmnij) $iP Sjp - 

g 47r V47T/ 
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Figure 1: The comparison of contributions Skin and Sf nt as the functions of the path along x axis. 
The instanton action is (3 = 18 and the source intensity is taken as e = g. The rotation matrix 
u ak is defined by the parameters n = (0, 1,0), <fi = n/2 (see Eq.(j2HJ)). Two upper curves (in the 
vicinity of coordinate origin) demonstrate the Skin behaviour whereas two lower ones demonstrate 
the S* nt behaviour. The solid curves correspond to the instanton and the dashed ones correspond to 
the anti-instanton. The source e is oriented along the third axis in the isotopic space. 



I g 3b k bl ^ e 3b k bl 
~. — Vk4m(Vlin Sj mn + fjlmn T n ij)£ a SiUJ a 5jLO — 5° fjkAmVlnc Zmcijn ^iUJ a 5jU + 

g 47T g 47r 

1 e 

i _a3c — — — /f. x^abdi bm c dn .bnr ,dm\ ^dran r ,adlr ,ck • 

gAir 

+ (y^j {Ski Uij - fj Mm fj Un V ijmn )(5iU ak 5jU al - 5iU 3k 5jU 31 ) + fj k4l V Uj 8iU 3k 5jp + 
1 e 



(17) 



a [{VlimVkmn Vk4mVlmn)^ijn + {VlArnVkin VkimVljn) -^mni "I" ^VlimVkin-^mnj]^ £ $jP + 

gAir 
1 e 

i _a3c ,cfc— — — o \^abd( .bmc .dn ,bn r ,dm\ _dmnr ,adi r „ , 

gAn 

+ -^~ VdkVdm E knijl [e 3ab (uj a %uj bd - uo a %uj bc ) - e ac %uj 3a ] 5 jP . 

Here we are studying an impact of fully developed Coulomb field on an instanton and that is why 
the above formula includes the spatial indecies only of the deformation fields S^p, 5 u uj. The overt 
but cumbersome representations of the tensor components are given in Appendix. In this paragraph 
we are interested in their asymptotic values at A — ► 00 only and it is convenient for further analysis 
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to use the dimensionless coordinates according to the following change Aj/p — > Aj . The terms of 
0((5p) 2 ) are presented by three tensors. One of those is the third rank tensor 



^ijk — $ijA k Tji + 5ikAj £ 2 + £jfcAj £ 3 + AjAjAfc E 4 

49vr 2 In A 

"a 



£1 ~ -2tt 2 , £ 2 ~ 2tt 2 , S 3 ~ — , £ 4 ~ 



where Aj = is the unit vector; another one is the entirely symmetric tensor of the second rank 

Qa = $ij Qi + Q2, (18 

with the components 



,lnA 

"a 1 

and eventually the entirely symmetric tensor of the fourth rank 



Q x ~ 6vr 2 — - , Q 2 ~ ; 



Rijki = {SijSki + SikSji + SuSjk) Ri + (SijAkAi + S ik AjAi + <^A,-A fe + 

. . . . (19) 

+ 5 M AiAj + SflAiAk + 5 jk AiAi) R 2 + AiAjA k A t R 3 , 



where 



R, ~ vr 2 ^y , R 2 ~ , i? 3 - . 



,lnA 
"A 1 

If the particular component of any tensor is small comparing to the leading asymptotic term and 
is out of discussion we imply zero value for it and don't show its asymptotic form. All the tensor- 
functions, excluding the components £1, £ 2 of Eyfc, are decreasing with the distance increasing. An 
appearence of these nonzero asymptotic values at A — > 00 results from the contribution of interference 
term originated by the product of source field strength G^B) and the G^A) component of second 
order in (Sp) 2 . Such behaviour is generated by the substitution where the parameter fixing above 
mentioned radius of multipole expansion L — ► 00. 

The coefficients of the 0((5u) 2 ) terms are constructed by two entirely symmetric tensors S and 
T of the third order. They are 

Sijk = (5ijA k + 5 ik Aj + 5 jk Aj) Si + AjA^A*. S 2 , (20) 

with the components 

- 2 lnA 



S 1 ~- — — T , S 2 ~0, 



7T~ 

"3 A 2 

and the tensor T is analogous to Eq.()20|) with similarly looking components 



vr 2 In A 

Besides, Eq. lj 17)1 contains the fifth rank tensor Z in the following form 

Zij k lm 

[Ai(5j k 5i m + Sji5 km + Sj m 5ki) + Aj(5ik5i m + SuSj m + Si m S k i) + 
+ A k (SijSi m + SuSjm + S im 5ji) + Ai(5ij6 km + S ik 6j m + S im Sj k )} Z\ + 
+ A m {bifi kl + 5 ik 5ji + 5 a 5 jk ) Z 2 + A m [8ijA k Ai + 5 ik AjAi + 5uA j A k + (21) 
+ 5 jfc AjA, + 5jiAiA k + SkiAiAj] Z 3 + [^AjAfcA, + 
+ (fjnAAfcA, + SkrrAiAjAi + ^AjA^Afe] Z 4 + AiAjA k AiA m Z 5 , 

their components have no singularities and are going to zero in the limit A — > 00. There exists the 
contribution of entirely symmetric tensor of the fifth rank $ in Eq. ()17j) which has the form analogous 
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to Z with the components $x — ^2? ^3 = I n the Appendix below two first components are 
denoted by $1, another two as $2 and, finally, the component $5 is denoted by $3. Unlikely the 
tensor Z the component of tensor $ having two time indecies ($4/34^) provides the terms of kinetic 
energy type with non-zero contribution. We demonstrate two nontrivial components of construction 
of the .Sijfc-type (Eq. (f2Uj) ) which are denoted as Qf\ in Appendix. As tensor Z the tensor $ does not 
develop the singularities in the limit A — > 00. The other terms of Eq. (jl7|) are given by the entirely 
symmetric tensor of the second rank U similar to Eq. (jl8|) with the components 



In A 

A 2 



U x ~ vr 2 — , U 2 ~ 



and the entirely symmetric tensor of the fourth rank V of the Eq. (jl9|) type with the components 



vr 2 In A 
TTx 2 



Vt--— — -, V 2 ~0, K 3 ~0 



In distinction with the dilatation term E in the contribution of kinetic energy type Eq. (j!7)) the term 
generated by the rotations in the isotopic space has no singular contributions (they compensate each 
other at the stage of intermediate calculations) and in output the tensor components are decreasing 
with the distance increasing. 

For the mixed contributions of 0(5p 5u) we have the tensor with the component approaching a 
constant value again 

^ ijk = (SijAk + 5 jk Ai) ^! + 6 ik Aj ^ 2 + AiAj-Afc * 3 , 
^1 ~ , ^2-0, #3 ~ 47r 2 ; 

(the singularity in the third component is related to the term jA of the action Eq.(Q)), and two other 
entirely symmetric tensors W and X of the third rank (of the Eq. ()20|) type again) with the following 
components 



Wt ~ -2vr 2 ^ , W 2 ~ , X x ~ -^-=^2- , A 2 ~ . 

The entirely symmetric tensor of the fifth rank 6 available in Eq.(jl7|) is constructed to be similar 
to $ with the nontrivial components 61,2,3, ®\\ = ®4ii34j- And, finally, the fifth rank tensor H takes 
the form of Eq.(|21|). Both tensors and 5 have no the singularities in the limit A — > 00. 

Now collecting the asymptotic values obtained all together we have for the additional contribution 
to the kinetic energy at A — > 00 (within the logarithmic precision) the following equation 

S? nt -> 2^%^ (—AjSu + AiSij) 5 lP 5 jP = , (22) 

which implies that any singularities available in the components of tensors do not produce any impact 
on the observables and, hence, the simplified procedure of the Ritz method with the parameter 
L — ► 00 is quite applicable. 

The complete analysis of the PP deformations including the contribution of the term S* nt Eq.(JT7J) 
may be cause for concern. However, the approximate procedure of searching the optimal deformations 

could be suggested just for the parameters interesting for applications in the IL model {(3 = ~ 

12 4- 18) [T3]. Let us address Fig.l where the contributions Skin and Sf nt for the single instanton are 
compared as functions of the path along the x-axis. We suggest for the instanton action (5 — 18 and 
the point-like souce intensity is taken as e = g. As to the rotation matrix it is defined by the vector 
n = (0, 1, 0) which fixes the PP rotation in the isotpic space to the angle ip in the following form 

u ab = n a n b + P ab cos ip - e abc n c sin p , (23) 



In A „. „ 4vr 2 lnA 
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where P a b = 5 a b — n a rib is the appropriate projection operator. In the example under discussion it is 
taken ip = ir/2. If one neglects the contribution of S* nt (comparing it to the contribution of Skin) the 
deformations 5p, 5u could be calculated in the true form (see Eq. (JH2j) below). It is clear from Fig.l 
the application of this approximation is well justified at small distances from the point-like source for 
both the instanton and the anti-instanton. The contribution S^ nt for the anti-instanton in the region 
of instanton size order is rather large formally but in this case the contribution of the term of kinetic 
energy type is small compared to Sf nt and the solution may be also used for the rough quantatative 
estimate. In what follows all the figures are obtained for the same parameters of rotation matrix and 
the same path as in Fig.l. 

i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 




■ ■ ■ i i i 



■ ■ 



A/p 



8 



10 



Figure 2: Simile of the contributions to Si nt . The path and parameters which define the rotation 
matrix are the same as in Fig. 1. Two upper curves show Sf t and two lower ones demonstrate 



Si nt . The solid curves correspond to an instanton whereas both dashed curves correspond to an 
anti-instanton. 



2.2 Analysing interaction of the undeformed instanton 

Here we are going to analyze the other interaction terms in Eq. (jl6|) . The second term there hinges 
directly on the asymptotic value of instanton size p and orientation u ak and describes the interaction 
of undeformed instanton with the point-like source 

St. = \ i *« + {£) J + {£) K U ,V , (24) 
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It is interesting to notice the first term with the coefficient 



I { = 8 tt 2 A { (D- 1 - A' 1 ' 

where D 2 = A 2 + l is entirely generated by the term jA of initial action Eq.(£Q). And the contributions 
of chromoelectric and chromomagnetic fields to Sf nt in two terms 

2 G%(B) G%,{A) + 2 G%{A) G%{A, B) , (25) 

after integrating are cancelled. If it occurs the instanton field is a weak perturbation of the Coulomb 
background (just the CDG approximation) then the first term becomes dominant and one can obtain 



\ J dx 2 G%{B) GIM) = £ — co 3k ^(2D 3 -2 A 3 - 3D + D' 1 ) 



Turning to the limit A — > oo we get the well known result for the dipole interaction of an instanton 
with an external field 

\jdx2 G%(B) G%(A) -> — f} Ui oo 3k Et , 

if the transition to the Minkowski space is performed and the dimensional variables are restored (see 
remarks after Eq.(jHJ). The vector Ei above presents the source field in the following form 

e Ai 

47T A 3 ' 

The function 

t-o [a 2 v 2 -y 2 

J ~ Z J dy (y 2 + l) 2 | ?/ + A| 2 ' 

and tensor 

K kl = 2 / dy 2 2 — ——2 , 

J (r + 1) 2 \y + A\ 2 

cannot be integrated in the elementary functions although numerically can be analysed comprehen- 
sively. Their asymptotic values at A — ► oo look simply 

5tt 2 1 

'J 



2 A 2 ' 

and for the components of the second rank tensor Eq. ()18jl as 



7rM 

2 A 2 



K x ~ —-2 , K 2 ~ . 



Now taking for the action of single instanton, for example, /? = 18 and for the source intensity e = g 
we find out the contribution of the terms proportional e 2 in Eq.(|24p at small distances (A < 1) is 
noticeably larger than the first term contribution. The behaviour of the curve corresponding to Sf nt 
in Fig. 2 allows us to make such a conclusion. If the contribution of the first dipole-like term in Eq. (j24j) 
is the leading one the dashed curve corresponding to the anti-instanton should run completely in the 
region of negative values. Further on in this paper we face another interesting result related to the 
slow decreasing of these functions while applied to IL. 
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2.3 Linear deformations 

Apparently, the interaction terms which allow us to detach the optimal deformations of PP are of 
major interest in this paper. Those terms are just linear in 5p and 5uj and could be given in the 
following form 

Si t = k n 5 nP + M* k 5 n u ak } (26) 

where 

An = - ^ fj k4i u 3k A in + (j-^j B n + (j-^j u 3k u 31 C knl . (27) 

The tensor components are cited in Appendix and here we give their asymptotic values at large 
distances. For example, for the second rank tensor A which is similar Eq. (jl8|) these values appear as 

A Y ~ -2vr 2 4 , A 2 ~ 16tt 2 4 • 
A A 

The vector-function B>i spanned on the unit vector A has the form 

* 28tt 2 In A 

B>i = Ai B , B ~ — 5- 

3 A 3 

and two components of the third rank tensor C similar Eq. (j2()j) are 

4vr 2 lnA „ 

Finally, the term in Eq. (j26j) related to the rotations in isotopic space is defined by the following 
tensor- functions 

1 e 1 e t 

M% = - —fjkAi D in ^ ~. £ a ° W° (Vk4mVlni E mi + fjk4jVljm Fmn) + 

g Air g 4tc 

1 e 

+ - — [e a3b (v bc (fj c4 jf}kii - VMjfjcu) - 5 a3 E kbc fj Mj f] ci i} H jtni + 

q 47T 

1 e (28) 

+ - — [(cA 3c - S 3a 5 bc )S kl - (u ab u 31 - S 3a S lb )S kc - e a3d u db e kd } fjuaVc^ms W + 

9 47T 

+ (^) 2 (co al - co 3l 5 a3 ) (OJ kl - P lkn ) . 

Then the corresponding asymptotic values when A — > 00 for two components of the D tensor similar 
to Eq.flHJ) are 

D 1 ~ , D 2 ~ 8vr 2 ^ . 
For the components of tensors E and F similar to Eq. ()2())l we have 

Ei ~ vr 2 l , £ 2 ~ , Fx ~ 2tt 2 1 , F 2 ~ . 

There is nothing special to get for the vector function 

0, = A,0, 0~-4tt 2 ^ 

and for the components of the third rank tensor P similar to Eq. (j20j) 

2vr 2 1 8tt 2 1 

Pi —~ , P2 ~ 



3 A ' z 3 A ' 
12 



The fourth rank tensor H reads 



Hijki = {Sij5 M + dikdji + 8 u 8 jk ) Hi + (<%A fc A z + S ik AjAi + 8 jh AiAi) H 2 + 
+ {5 a kjk k + 8 M AiAj + SflAiAk) H 3 + A^-A^A, H 4 , 

with the following asymptotics of its components 

H - --— - --— H -> i— -> -- — 

1 36 A 3 ' 2 12 A 3 ' 3 2 A 3 ' 4 2 A 3 ' 

For the space components of entirely symmetric tensor of the fourth rank Y similar to Eq.(|19|) we 
receive 

lvr 2 

Besides, the contribution of tensor with two time indecies Y±ij± (similar to Eq. (jl8jl ) is non-zero as 
well. Its components Y\ \ in the limit A — ► oo are the following 

y44 _^ \± y44 ^ 33^ 
1 6 A ' 2 16 A ' 

Now making use of the definitions of A n (EqEZTj) and M® k (EqEHJ) we find the interaction Sf nt con- 
tributing the terms of the second order in the parameter to the action. Those terms are quite 

comparable with the contribution of Sf nt in the regime interesting for physics when the source in- 
tensity e is of the same order as g (the factor 4tc is obviously compensated while integrating). The 
suppression of multipole contributions is expected to come from the geometric factor which emerges 
while integrating over the azimuthal angle and for the field configurations rather homogeneous in 
space, approximately it behaves like 2k\- 1 wnere ^ ^ s ^ ne multipole type. Fig.2 where the contri- 
butions of Sf nt (two upper curves) and Sf nt (two lower curves) are exhibited manifestly demonstrates 
the role of interactions of deformable PPs (solid curves correspond to the instanton and the dashed 
ones correspond to the anti- instanton). 

3 Developing an approximate scheme of deformations 

Now we are well equipped to find out the optimal instanton deformations. As mentioned above the 
corresponding formulae can be obtained analytically if we neglect the contribution of S^ nt Eq. ([T7Jl 
which is distance dependent. Optimizing the quadratic form in the action Eq.([15|) we receive the 
following system of algebraic equations 

k 5 iP + v* 5ju al + Ai = , 

trg 5jU al + vg S jP + v mij e abc v 2 (2 ui bm 8^ ck - u b %u cm ) - (29) 
-v 2 fjkije^uj^uj™ - 2«i e kmn e abc uj bm 5 % uj cn + \v 2 uj ak uj cn b { uj cn + Mf = . 

In order to resolve it we introduce a new variable 

y jki = u a j5iU ak _ 

Recalling the orthogonality property of the rotation matrices and invariance of the group measure 
Eq.(|T3J) we have 

<7i u al y m + a 2 5 ki u al yi nn - o 2 uj al y lik + 2v 2 e mij e dml u ad y lkj - v 2 e mij e dkl u ad y lmj - 

„ .ad~ Itij o„, „ „ .Oid„ Ini , ,ak„ nni , ~ak r „ , i\^ak n 

-v 2 ekije dn i lo y J - 2v 1 e kmn e dm i u y + -v 2 u y + v {j djp + M i = . 
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Multiplying it by the u matrix and calculating the products of the antisymmetric tensors we discover 
the equations 

ox f jk - 02 y ikj + 2 Vl y^ k + v 2 y* ki + v 2 y k ^ + 

+S jk (a 2 y mn - v 2 y nin + v 2 y nm ) - v 2 5 kl (y^ nn + y nnj ) + (30) 
+5 jt 



-v 2 y nkn + v 2 y knn + (-v 2 - 2 Vl \ y nnk + vg 5 n p + uj m M a k j = . 



The solution of this equation system is essentially based on the presence of the convolution of inde- 
pendent variable y^. It is clear these equations are not contradictory if the following constraints on 
three possible convolutions y %kk ) y k%k ^ y kkl are obeyed 

(a, + 2a 2 ) y lkk + (2v, - 3v 2 ) y ktk + (^-v 2 - 2v^j y kki + u ai M£ k + vf n 5 n p = , 

(a 2 + 3v 2 ) y lkk - ((T 2 + 3v 2 ) y kik + (a, - Av x + Av 2 ) y kkl + u ak Mf + v kk 5 n p = , (31) 

(a 2 + 2v x - v 2 ) y lkk + (^ - 2v 2 ) y kik + (| - a 2 - 2v^j y kki + cu ak M^ + 5 n p = , 

For what follows it is practical to analyze a more general situation of the system with arbitrary 
constraints. The first five terms of Eq.()30p which do not contain the convolutions might be presented 
using the following auxiliary tensor 

s ipP,kqi _ Sl Sjk5p q 5j3ry + s 2 5jk5jry5f3 g + s 3 5j q S P k5f3 y + s 4 5j g 5 pj 5pk + S5 5j^5 P kSp q , 

with the components Si = a±, s 2 = —a 2 , S3 = 2vi, S4 = S5 = v 2 . The inverse tensor has the following 
form 

T ioa;jpf3 _ ^ §..§ op § af3 + T2 dijSofidap + T 3 5 ip 5 oj 5 af3 + T 4 5 ip 5 p5 a j + T 5 5 il3 5 oj 5 al 3 + Tq 5^6^,6^ , 

and it is valid for this tensor r %OOL ^ v " s^ v ^ kqi = 5ikS oq 5 ay and its components are 

11-199 25-31 35-37 



r i = — f; — > r 2 = f; , ^"3 



D D D 

109 7 • 23 

r 4 = r 5 = — , r 6 = — , D = 13 • 83 p . 

In principle these 'strength' coefficients could be presented in the terms of corresponding combinations 
of Oi and V{. However, they look like too much sophisticated and it is the main reason for showing 
their overt expressions here for PP in the singular gauge. It is clear that the relations obtained are 
enough to get the solution of the system Eq- lj^Uj) . The convolutions are extracted from the constraint 
equations Eq. pTf 



yk 

V y k 





136 


24 


8 


ikk \ 

kik 1 

kki J 


25 p 
16 


175 p 
24 

ll/ 


5 p 
8 




25 p 


175 p 


5 p 



u ai M ak + ^ _ 2v 2 )5iP 

u ak Mf + 3 Vl 5ip 
V oj ak M^ + («i + 2v 2 )5 iP 



When it is used in the equations for dp we obtain 

S3 28 19 2 

" 200 P Kp ~ 25 + 25 ^ Mfc ~ 5 + K = ° " (32) 

Making use of the tensor r which is inverse to the tensor s of the system Eq. (j30|) we find out the 
solution y and reconstruct the solution 5iU ak = io a iyi k% . The overt formulae are too cumbersome and 
are not shown here. 
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Figure 3: The norm of instanton size deviations \Sp\ for version with the parameters mentioned in 
Fig.l caption. The dashed line corresponds to the anti- instanton. 

4 Analysing numerical results 

The most interesting question we are going to answer here is what are the distances where the 
application of the perturbative approach developed is still valid and reliably justified. Actually, the 
formal criterion could be grounded on the deformation smallness on the scale of characteristic PP 
size which corresponds to the distance of order \y\ ~ 1 in the dimensionless variables. Fig. 3 illustrates 
the behaviour of the norm of instanton size deviations 

\Sp\ = max \5ip\ , i = 1, 2, 3 , 

i 

as the function of the path. That is clearly seen the perturbative approach is completely relevant 
at small distances and distances commensurate with several instanton sizes. Besides, we could also 
conclude the deviations from asymptotic PP size are strongly dependent on the mutual positions of 
PP and point-like source in the isotopic space. The example of Fig. 3 is one of the most rigorous 
limitations in \Sp\. If one orientates an instanton along the third axis of isotopic space collinearly 
with the point-like source the deformations get considerably smaller. This configuration in isotopic 
space is characterized by the prevailing contribution of rotational mode to Sf nt Eq. f|TH|) comparing to 
the dilatational one. However, if one orientates PP along the direction of point-like source vector, 
the corresponding contribution of dilatational mode becomes prevailing, i.e. it is more advantageous 
for PP to rotate in the isotopic space being correlated with the point-source. Fig. 4 displays the total 
increase of instanton action (S — (3) (solid curve) and anti-instanton (dashed curve) as function of 
the path. The upper solid line shows the contribution of Sf nt from Eq. fjlfijl in the proximity of the 
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coordinate origin. Let us notice here that at large distances both the instanton trajectory and the 
anti-instanton one are running above Sf nt which is just the action corresponding to non-deformed 
instanton. It means the crumpled PP develops the action maximum in this region and, hence, the 
paths in the functional space might exist and the PP action might be decreasing while going along 
them. 

In order to appoint those we consider more general problem if determine the extremal value of 
the functional Eq. (jl5|) with nine arbitrary constraints 

y ikk = cii , y kik = h , y kkl = q . 

The action Eq. (jl5|) with the Lagrange multipliers A, p,, v might be transformed into the form 

S^S = S + - a,) + ^{y kik - b t ) + v { {y kkl - a) . (33) 

In fact, it is not easy to diagonalize directly the quadratic form for this system. Nevertheless, the 
useful information on the general construction of this function can be extracted. Trying to do that 
we rewrite the part related to the rotations in the isotopic space 

a i P a;j P p = °_i s..$ pq s a/3 - ?1 5ij8 pl3 5 aq + v 2 5 iq 5 p p8 aj + v x 5 iq 5 pj 5 al3 + (34) 
using the multiplet expansion in the form 

yjqP _ y{jqP} + Jjqfi + y{j[q}/3] + y ij{q]ft 

where the combinations symmetric in indecies are denoted by the curly brackets and antisymmetric 
ones are given by the square brackets. Then the terms without convolutions can be written in the 
diagonal form as 

S y = y ipa 



ri y {ipa} + r2 y [ipa] + r3 y {i[p}a] + r ^ y [i{p]a} + 



(35) 



where the eigenvalues are defined as 



o-iG-2, 13 a x <T 2 1 

n = Y - Y + v 2 + Vl = -p, r 2 = Y + - + v 2 - Vl = -f3, 

ai 02 11 a x a 2 1 

r 3 = h v i — v o = — p , r 4 = v 2 — vi = -p . 

2 2 1 2 24 2 2 4 



(36) 



As all Ti are positive we find every separate multiplet to be stable and clearly distinguishable because 
of their contribution to the term of kinetic energy type. The system of equations with arbitrary 
constraints imposed looks like 

ttSip + v x y kki + v 2 y kik + v 2 y ikk + K = , (37) 
oxy ldk - a 2 y ikj + 2v iy > lk + v 2 y^ + v 2 y kij + 5 jk (a 2 y mn - v 2 y nin + v 2 y nni ) - v 2 5 ki (y^ n + y nn >) + 



+5. 
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-r>y U " +r.J>" , + (^ 2 - 2/-, 



+ Xi5 jk + fij5 ik + u k 5ij + v l k J n 5 n p + u ai M^ J = . 



Then the conditions for this system not to be the contradictory one might be taken as the equations 
for the Lagrange multipliers 

/13 \ 

3A + pL + v + {a x + 2a 2 )a + (2«i - 3v 2 )b + i—v 2 - 2vA c + c^Mf + {v x - 2v 2 )5p = , 
A + 3/x + v + (a 2 + 3v 2 )a - (a 2 + 3v 2 )b + (a x - A Vl + Av 2 )c + uj ak hlf + ( Vl + 2v 2 )5p = , 
A + p + 3v + {a 2 + 2v x - v 2 )a + (a x - 2v 2 )b + (j - a 2 - 2^) c + u ak M^ + 3 Vl 5p = . 

16 



The equation for the size deformation reads 



K5p — v 2 a + v 2 b — v\c + A = . 

These equations provide us with A, /x, v and putting those in Eq.(pT7j) we come (using the inverse 
tensor r) to the solution for the system with constraints. At large distances from the charge where 
the interaction is negligible or in the particular case of the interaction absent (the deformed PP only 
is considered) the Lagrange multipliers A, n, v and solution y^ h are the functions of constraints a, 
b, c only 




83 

I 

"a 



43 
353 
571 



49 

"377T 
571 
3 D L 

557 



3 D[ 3Di 




where D x = 5 ' 9 ^ 32 . 



ijk 



+ 5 Jk (A n Oi + A l2 hi + A 13 a) + 5 lk (A 21 a, + A 22 b, + A 23 Cj ) + 
+ 5 ij {A 31 a k + A 32 b k + A 33 Cjfc ) = , 



moreover, the matrix A is just a projector necessary to excrete from the solution y nine constraints 
imposed A n = A 22 = A 33 = A 12 = A 13 = A 21 = A 23 = A 31 = A 32 
action Eq.()33J) becomes the quadratic form of three constraint vectors a, b, c i.e. 



yg. Then as the result the 



S = — — (Aa + /j,b + vc) . 

Two eigenvalues of this form are positive and one is negative. The direction where going along the 
action is getting smaller is given by the following equation 



a 



0, 3-571 6- 557 c = 



(let us remind once more that these result was obtained considering the instanton in singular gauge, 
and we omitted their representations while the corresponding components of tensors cr,, Vi because 
of their very cumbersome forms). We did not manage to diogonalyze equations (|30|) manifestly 
nevertheless imposing arbitrary constaints we demonstrate the part of eigenvalues of quadratic form 
is negative and PP can manifest itself as a 'sphaleron'. In principle studying the deformed PPs out 
of the perturbation theory developed here could allow to understand transition of nonperturbative 
field configurations into the perturbative ones. Usually such a transition is associated with instanton 
anti-instanton annihilation only and justifies the criticism of simple superposition ansatz in the IL 
model [0]. From this view point we demonstrate that already in one particle sector there exist 
the ways in the functional space which inter-relate the solutions from different topological classes. 
However, it is a well known fact that the solution with nontrivial topology in spite of the loss of 
action is characterised by large statistical weight. 

Dealing with one PP we could study the deformations which leave the solution within the same 
topological class if the condition of topological charge conservation SN = is imposed. In the 
approximation of quadratic deviations we receive the following result 



3 A A 1 c 



ak 



ak 



(38) 



where G a 



1 fia 



It leads to the conclusion that the dipole-like deformations (the expansion 



Eq-0 is meant) are forbidden if the toplogical charge is not changed. 
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As an example of the system with the constraints imposed we consider the substitution with four 
non-zero components (which is dictated by the particular form of matrix uo) 



Sip = x , Situ 13 = 5iuj 22 = y , Situ 31 = z . 

Besides, the contribution of terms Sf nt could be analysed in this example. The solution is determined 
by minimizing the following quadratic form 

S-(3= — x 2 + ^y 2 + ^-(3 z 2 + Fyz + Cxy + Dxz + A lX + Ey + Gz , 
2 2 lb 

where 

1 = f T ^ SlU+ (i) (2^-^211-^3311), 

f = ^/5±^£ (2T m - 5 133 - 5 122 )+ (£) 2 (2f/ n - ^211-^3311) ± 

± ~ 2 (Z 1 -Z 2 -Z 3 + Z 4 ) , 
#4tt 

C = l-(3 ± 2 (2W 1U -3X 221 -X 331 )±-^- 4(S 1 -S 2 -S 3 + S 4 ), 

12 g An g Air 

D = -^-/3t-^*iii±-^2(H 1 -H 2 -H3 + H 4 )± 
24 g 47r g An 

± 2 (39f + 9^ + 219! + 129 2 + e 3 ) , 

5(47T 

13 1 e 

F = — /3±-— 2 (3$f + $^ 4 + 21$ 1 + 12$ 2 + $ 3 ) 
24 gA-n v 1 2 y 

E = M x 13 + Ml 2 , G = Mf , 

here the upper sign corresponds to the instanton. The instanton trajectory is depicted by the dotted 
line in Fig. 4. 

The topological charge change for the trajectories under consideration is plotted in Fig. 5 and 
allows one to conclude that the interaction does not take the solution out of the topological class 
with the charge \N\ = 1. The quantity defined by Eq.fjHBJ) characterizes the deformation extent of 
aproximate instanton solution only, also one should pay interest to topological charge of superposition 
solution on the whole. 

Another interesting remark as to the possible constraint conditions comes from their optimization, 
i.e. from the analysis of nine constraints obtained for the action extremum available. It turns out 
the initial trajectories of Eq. (}3T|) are just the case and the Lagrange multipliers for this extremal 
action are trivial (A = fi = v = 0). 



5 Analyzing the IL model 

The results obtained above allow one to conclude the perturbative approach developed is applicable 
to effectively explore the external field impact on PP no matter how large (or small) the distances 
are. One should expect the substantial suppression of instantons at the distances of average instanton 
size order in IL but the Coulomb-like interaction extends to the distances much larger, of several 
characteristic instanton sizes. 

Apparently, the configurations of crumpled instantons may come about very essential for the IL 
model. The entropy of crumpled pseudoparticles and their contribution to the functional integral 
are of the same order as for the habitual instantons. (The entropy of a field configuration is implied 
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Figure 4: The increase of action (S — (3) for the crumpled instanton and anti-instanton (dashed line). 
The parameters of the rotation matrix and the path are mentioned in the caption of Fig.l. Two solid 
lines correspond to the instanton and the upper one (around the coordinate origin) shows Sf nt . The 
dotted line shows the instanton trajectory with the constraints imposed. 

as the log of the volume in functional space occupied by similar configurations.) An importance 
of searching the similar trajectories which have the large entropy (even losing in the action) at 
calculating the functional integral is well known (see, for example, (El). In fact, the condition 
5S = for searching the optimal crumpled configurations operates with the interaction contributions 
coming from exponential factor only and, as can be shown, is a valley equation in quasi-clasical 
approximation. In order to pretend to higher precision and to take into account so-called quantum 
interactions one should deal with the valley method (see for example ^3] where the similar practical 
method is given). However, the discussed example of dipole-like deformations Eq.(7) demonstrates 
obviously how sophisticated this configuration can happen to be even in this comparatively simple 
case. 

Another remark worth mentioning here concerns the behaviour of the pseudo-particle deforma- 
tions at large distances. Studying the asymptotic behaviour of interaction terms one easily finds out 
their Coulomb-like character, i.e. even in far distant area the pseudo-particle field may be noticeably 
polarized. In principle, if it turns out possible in some way to detach the long-range component on the 
background of inter-instanton interactions (for example, for dilute gas of very cooled (anti-)instanton 
configurations) the chance may appear to measure the coupling constant g on the lattice directly if 
an external field of intensity e (which is controlled) is applied to its boundary. The dilatational and 
rotational deformations may contribute to the observables only as the corresponding averages (5 n p), 
{Suf b ) which are defined by two independent parameters e 2 and §. It looks feasible to construct such 
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a linear combination of those two averages which depends on the intensity e only and certainly should 
be under control during the measuring process whereas another could provide us with the magnitude 
of g. The great experience collected in lattice QCD to work with the boundary conditions imposed 
[TH] gives us the hope to succeed on this way applying the cooling technics [THj which isolates the 
instanton configurations successfully and incidently they look quite deformed [T7j . 



5.1 Crumpling instantons in the IL model 

Here we are trying to suggest a fairly simple approach of taking into account the crumpled pseu- 
doparticle configurations in the functional integral. If we pretend to make this estimate calculating 
precisely the exponential factor of generating functional only then their contribution could be found 
in quasi-classical approximation just in the same way as it is done for the undeformed pseudo-particles 
(see the designations in [T2] ) 

1 N . j 

z = e w n J v^ Cnc 13 e ~ s ' dii = dZi dpi dujt - (39) 



N ■ i=l 



As supposed in the IL model this integral is saturated by the configurations of characteristic size 
p. In the approach under consideration S means an action including the deformations in the above 
formula. Thus, adding the contibution of the crumpled pseudo-particles of the similar characteristic 
scale we are able to work within the same precision studying the IL reaction on the external impact. 
The variational Ritz method as defined by Eq.© looks to be pretty hard and incomplete from the 
viewpoint of the integration in the functional space. We have permitted the velocities and higher 
coefficients of multipole expansion only to vary there but not the major expansion term that is the 
average instanton size. In view of the functional integration expanding the average pseudo-particle 
size varied seems more natural. Then we have for this particular expansion the following expression 

r> i x / x . dp{z) d 2 p{z) y^ y v 

J 2ta (s > *)=p(*) + — y, + — — + ..., \y\ < L . (40) 

and the similar one for Q. The generating functional becomes now 

I SL f dz,; duji 

N 



J D[p(z)) D[u{z)\ E ^[ / C Nc P™< e~ s . (41) 



Remembering that the IL approach, as formulated, takes into account interactions between the 
pseudo-particles we should expect in Eq.® the pseudo-particle as the field B, too. In order to 
simplify the calculations let us suppose for the moment the (anti-)instantons undergo the insignificant 
changes of size because of an impact of external source on IL, i.e. p(z) = p+5p(z) and colour instanton 
orientation gets u>i(z) = u t + 5ou(z) and for anti-instanton it is Ui(z) = Ui + Su(z). The first terms of 
uj correspond to the anticipated fixed directions in colour space. The field strength of single instanton 
looks like 

G; U (A) = -~ g ^\z)M, a M v , fj ka p {y2 f 

with uo changed for uj for the anti-instanton. When we have the instanton-instanton 'molecule' the 
mixed component of field strength G(A, B) reads 

4 „abc .fen/ \ ,cfc/ N t=i * — — \ yi*y2f p 2 (^l) p 2 (z 2 ) 



G a ^(A,B) = - e abc oj bn (z 1 )uj ck (z 2 ) (r] n ^r] kl/a - r] nvi ^ km! 2 2 



where yi = x — Z{. It is clear all the combinations of colour matrices for the instanton-instanton 
interaction owing to the property Eq. (fT4*|) give the unit matrix if we limit ourselves with the contact 
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type of interaction only. This approximation seems well justified since the integrals of our interest are 
quickly decreasing for the far distant pseudo-particles \z\ — z 2 \ ^> p. For the instanton-anti-instanton 
molecule such a component is given by 



G%,(A,B) = - e abC U bn (z 1 )uJ Ck (z2) (Vnfj.jVkua - VknaVnuyj 2 2 



Vlvl (y 2 1 +p 2 (z 1 ))(y 2 + p 2 (z 2 )) ' 

The properties of the 't Hooft symbol make it possible to demonstrate that the product of matrices 
uj and uj gives the unit one in this case as well. Then the averaged integral for ineracting pseudo- 
particles may be written in the following form (see Ref. for designations) 

(S(AB)) = j^U int (A,B), 

U int (A, B) = I** ( 7 ^-^) 8 )^)^) ^ 



g 2 N 2 - 1 J y\{y\ + p\z x )f jfi(y* + p\z 2 )f ' 

And integrating over the positions of single detached pseudo-particle one may receive an estimate in 
the form of contact term 

/ dzi u int (AB) ^-^e P 2 {z 2 )p\z 2 ) ^^-e p 2 p\z 2 ) , 

J 9 9 

where £ 2 = ^^~luT^~T- (Obviously the interaction integral is independent of the colour orientation 

c 

because of the IL isotropy in colour space.) Besides, we should take into account the term of kinetic 
energy type Eq. fjllj) for each pseudo-particle. As the result, it has been shown in jjj, the generating 
functional (jH)) for such a sort of crumpled pseudo-particles takes the form of the effective Lagrangian 

for the scalar field p(z) with the mass gap M 2 = - ^~^ c q — — The dependence on the 

O K p 

rotational component (5u, Sou) occurs to be trivial. The estimate for kinetic coefficient k ~ 1.5/3 — 
6/3 has been obtained before and depends on the ansatz for the saturating configuration. For example, 
for pure gluodynamics with the parameters iV c = 3, f3 ~ 17.5, pA ~ 0.37 at k = 4/3 it is M ~ 1.21A. 
In this case we have the soft mass scale with characteristic for IL A -280 MeV. The result k = 0.9/3 
we have received in this paper corresponds to essentially harder mass scale, above 1 GeV. The IL 
excitation discussed in its quantum numbers could be identified as a glueball. 



5.2 Average source energy in IL 

In this paper we are studying the way in which an external field acts upon the pseudo-particle. 
If we don't fix the features of field B°(x) in advance the feedback of IL could be described self- 
consistently. Then the multipole expansion we used allows us to factorize the problem in a sense 
and construct the closed equation for the Fourier component of the field B°(p). That equation looks 
rather promising to calculate the Debay screening radius that is one of the most important physical 
parameters characterizing IL. However, even the description we have already elaborated for the 
(anti-)instanton behaviour in external field is pretty rich to provide us with interesting information 
on the IL properties and possibility to explore the various observables. We mention here one result 
only calculating an average energy of Euclidean source inserted into IL. As a relevant saturating 
configuration in the IL model one takes the superposition of (anti-)instanton fields with the source 
field _B"added as 

N 

A«(x) = B«(x) + 5»; 7i ) , (43) 
i=l 

where 7, = (pi,Zi,u>i) denotes the parameters describing the i-th instanton. It is seen in Fig. 4 that 
at the distances larger than 2A/p the IL density becomes practically equal to its asymptotic value 
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-0.0 




Figure 5: Behaviour of the topological charge 5N for the crumpled instanton and anti-instanton 
(dashed line). The parameters of rotation matrix and path are the same as in Fig.l. The dotted line 
corresponds the trajectory for instanton with the constraints imposed, see the text. 

n(A) ~ n e /3_,s ~ uq because the action of each pseudo-particle there approximately coincides with 
j3, S — (3 > Sf nt . The quantity here we are interested in should be given by averaging S over the 
pseudo-particle positions and their colour orientations (taking all the pseudo-particles of the same 
size for simplicity) as 

<s> = n/^ /** s= 

-£f*+">W£(=)'('+¥). 

where L is a formal upper integration limit, V = L 3 X^ defines the IL volume, X4 is an upper bound 
of the 'time' integration, N denotes the total number of pseudo-particles and a is a source size value 
(on strong interaction scale, of course). Taking the asymptotic values of J and K functions and 
returning to the dimensional variables for the moment one instanton contribution to the average 
action can be written down in the following form 

where p is the mean size of (anti-)instanton in IL. The result is given in the form with the common 
factor X4 extracted because our concern now is the (anti-)instanton behaviour in the source field 
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background where the field is steady and the solution possesses the scaling property at any x 4 slice. 
In the limit N, V — > oo and at the IL density n = N/V fixed this result becomes 

(S)^EX 4 , E=^- + nf3L 3 + n^^p 2 L. 

An a p g 

Formally, the last term of E looks like a small correction to the gluon condensate (the second term). 
However, this contribution linearly increasing with L and proportional to e 2 has a different physical 
meaning of an additional contribution to the source self-energy 

_ _ 67r 3 e 2 n 

E ~ a L , a = n — — 5- p . 

P 9 

The tension value a for the IL characteristic parameters ~ where where .R is the mean separation 

of pseudo-particles, n = R~ A , (3 ~ 12, p ~ 1 Ge^ -1 [H],[I] comes about a ~ 0.6 GeV/fm (if one takes 
for the brief estimate for the source intensity e ~ g). This result looks fully relevant in view of the 
qualitative character of estimates which IL is able to provide. Moreover, such a value is in reasonable 
agreement with the estimates extracted from the potential models for heavy quarkonia, for example. 
If one intends to explore the magnitudes like (Si nt e~ Sint ) (which could model an effect of suppressing 
the pseudo-particle contribution in the source vicinity) in numerical calculations it becomes evident 
the linearly increasing behaviour starts to form at A/p ~ 3 — 4. An interesting phenomenon of energy 
increase with a distance increasing which we discovered above looks promising for the applications 
in IL. Indeed, this result teaches the source mass (we should treat an additional contribution found 
out just in this way) is unboundedly increasing (if the screening mechanism is not included) what 
could be interpreted as preventing the possibility to immerse a bare colour charge into IL. Our 
estimate of asymptotic energy of the Euclidean source put into IL shows the major contribution to 
the generating functional in quasi- classical approximation while all coupling constants are fixed at 
the scale of average instanton size p. 

Further we are going to exploit this interesting feature of IL for analysis of the colour dipole 
behaviour dealing with the point-like sources ed a3 , —ed a3 which are oriented along the third axis of 
isotopic space and placed in the points z% and z 2 , respectively. Then the strength tensor component 
of our interest looks as follows 



GUA, B) = 2f e a3c Q ck fj kia % . 2 p2 2 . ( - 




4?r ~ — y 2 {y 2 + p 2 ) y\ y + Al \ \y + A 2 \J ' 

where A 12 = z— Z\ i. The contribution to the average energy E at large distances which is additional 
to the Coulomb one is defined by the integral similar to that for the configuration with single source 
(the same coefficient but distinguishing average over the instanon positions), i.e. 

I d ~- 

When the distance between the sources I = |Ai — A 2 | is going to zero the field disappears and the 
final result should be zero. The integral may be dependent on two parameters L and I only. The 
dimensional analysis reveals the integral is a linear function of both parameters but the /-dependence 
only obeys the requirement of integral petering out at I — > 0. In order to determine the coefficient 
we find easily 

(here the contributions of three integrals are shown separately), where I = |Ai — A 2 | is the distance 
between charges. And finally we have for the additional contribution to the dipole average energy 
while in IL 

E ~al . 
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The non-singlet in colour states are strongly suppressed in IL by the factor of order of the screening 
radius ~ a R D . 

Finally, we would like to make one general remark which is not directly related to IL but seems 
to be of heuristical importance. We encounter the crumpled configurations of topological nontrivial 
solutions (pseudo-particles) not only while studying the impact of external fields on them or analyzing 
their interactions. We discover such configurations, for example, studying the simplest kinks in the 
problem of double well potential V(x) = —X(x 2 — rj 2 ) 2 (we are dealing here with an analysis in the 
Euclidean space). Let us consider the kink at finite interval (— T, T) where T ^> uu~ l if the latter is a 
characteristic kink size (uj 2 = 8Xr] 2 ). The approximate solution is given by Xk(t) ~ x as (t) = r] th^ 
and at the turning points x k (T) = :r fc (— T) ~ r\ the particle velocity disappears Xk{±T) = 0. The 
trajectory choice is dictated by the requirement that a particle should pass through the coordinate 
origin at zero time moment (the coordinate of kink center t c — 0). Apart from this (unique) symmetric 
kink the trajectories with the biased center which crosses the axis x = not only at zero time moment 
are of special interest. In the vicinity of kink center such solutions are approximately described by 

x c (t) = r] th ^^ 2 ■ Then, in order to construct these solutions we consider the time interval larger 
than the initial one (— T 2 ,T 2 ) where T 2 > T. There is the trajectory with the kink center located 
at the coordinate origin which we designate as Xk(t), i.e. x^it) ~ x as (t), and this particle starts 
and finishes again with the zero velocity x k (±T 2 ) = 0. Calculating the energy of this particle we 
find the magnitude which is slightly larger than for the trajectory Xk{t) and the particle climbing 
the saddle point reaches the points higher than the turning points Xk{±T) of the initial trajectory. 
Now, let us shift the time interval origin for the trajectory Xk(t). Thus, we receive the solution 
x c (t) = x k (t + T — T 2 ) with the baised center. For this configuration the particle starts to fall 
down at the time moment — T to the well with the zero initial velocity and crosses over the axis 
x = at the time moment T 2 — T. Then climbing the top again it does not stop at the time 
moment T and continues to move with a certain velocity. Hence, we need to install a hard barrier 
to reflect it compelling to move in the opposite direction. Strictly speaking such trajectories are not 
the equation solutions for the initial potential. The trajectory of particle finishing 'properly' could 
be constructed by a similar way, however, in distinction from the kink x^{t) this trajectory x c {t) is 
asymmetric. Its right hand wing does not copy absolutely the left hand one. Truely, the difference is 
quite inconsiderable (exponentially small) if t c <C T but becomes very essential when the shift of kink 
center is of the same order as the interval magnitude T. Generally, the particle could have non-zero 
velocity both at a start and at a finish. If we remember now that calculating the functional integral 
we are limited by the special requirement for trajectories to deal with only those which have both 
ends fixed. Our discussion above teaches it could be satisfied by breaking the equation of motion. 

Thus, we are forced to examine the inconventional kinks which are just the particular set of 
trajectories with the measure in the functional space proportional to the magnitude of interval 
/ dt c — 2 T. The kink-anti-kink configurations form a much more sophisticated set of trajectories 
(with complicated interaction between fragments) which are transformed into standard oscillator 
(perturbative) trajectories, with oscillation frequency Q 2 = 4 A rj 2 , while approaching the well bottom. 
Such trajectory should not be identified with the kink-anti-kink configurations. The instructive 
message of this example, at least, in the context of previous study is that the exact definition of a 
kink is of less importance compared to the necessity of taking into account properly the contribution 
of adequate trajectories in the functional integral. We mean the trajectories of considerable power in 
the functional space, for example, integrating over their deformations if those contribute significantly 
to an effective action. 
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Conclusion 

Since the instanton discovery a suggestive way to access an information on the nonperturbative QCD 
dynamics is to investigate the role of this coherent gluon fields as major objects populating the QCD 
vacuum. The development of the IL model allowed to obtain the estimates, the bulk instanton 
properties in vacuum and the lattice simulations have later supported these. Hovewer, the study of 
instanton behaviour in the external fields, as was understood, is of great importance for resolving 
the principle QCD problems but the results here are still rather hazy. It was a main motivation to 
have in this paper as a focal point analysis of the PP behaviour in the field of point-like source. In 
order to deal with such a task we developed the perturbation approach based on the variation of 
the instanton parameters. We have formulated the variational problem of searching the optimal PP 
deformations which is algebraically resolved by the Ritz method within the multipole expansion. The 
regions of parameter changes interesting for the applications in the IL model were especially analyzed 
in detail. We demonstrate the solution of algebraic problem might be essentially simplified and the 
overt formula for the deformation fields in dipole approximation might be obtained if one neglects the 
contribution to the term of kinetic energy type S K which depends on the distance to the perturbation 
source. The coefficients of the quadratic form in this term of kinetic energy type S K are related to PP 
itself and charaterize, in a sense, its 'pliablity' for changing the size and orientation in the 'isotopic' 
space. The terms responsible for the interactions between undeformed PP and PP with unchanged 
parameters with fixed point-like Euclidean source of colour field were calculated. In addition, it was 
realized the dipole-like contributions to Sf nt which are proportional to e/g occur noticeably less in the 
parameter region characteristic for the IL model at smaller distances than the terms of e 2 -order. The 
effective and rather simple method to take into account the 'crumpled' configurations at calculating 
approximately the functional integral. The description of instanton behaviour in the field of point-like 
Euclidean colour source developed in the paper clearly demonstrates that the valley topography in the 
functional space might be so sophisticated that makes it impossible to investigate the problem in any 
general formulation. Finally, in the framework of superposition ansatz the estimate of average energy 
of non-abelian dipole in the IL medium has been found out. This energy for the dipole in colour 
singlet state escalates linearly with the separation increasing for the point-like sources unscreened. 
The corresponding coefficient of this dependence develops the magnitude similar numerically to one 
as inferred from the lattice OCD calculations. 

Our work was partially supported by the Grants STCU #P015c, CERN-INTAS 2000-349, NATO 
2000-PST.CLG 977482. 



Appendix 

Here we present the integrable terms. We show only the components of the corresponding tensors 
which are mentioned in the paper The interacting terms linear in 5p are the following 
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The components of kinetic energy tensors which are the cross terms between 5p and 5u are defined 
as follows 



^2 = 
^3 = 



9 9 



A 3 
4tt^ 
A 3 
4tt 2 



-A 

9 



8 
9 J 



ln(A + D) 



, w - 28, v , 35 D /4A 3 
7 -L> 7 H L> 5 D 3 H h 1 



w 2 



A 3 

i6tt 2 r 2 



9 9 
8A 7 + 8D 7 - 28D 5 + 35L> 3 



6 ' V 3 
31 



-D 



2A 
15 



D 2A 



ln(A + D) 



A 3 
16tt 2 
A 3 



25 



(jD 9 _ A 9 } _ 1^7 ^_ D 5 _ 21 ^3 2 _ 

v ; 25 100 40 5 



A 



9 



8 40A 



ln(A + D) 



8 /a q ^ 36 - 63 . 21 „ 33^ 3 1 9 , /A 

- A 9 - L> 9 + — L> 7 L> 5 + —D 3 D + ^ + — n A + £>) 

5 V ; 5 5 2 8 4D 4.D 3 8A 



X 1 
X 2 



87T 2 

A 3 

87T 2 

A 3 



D 7 A 7 7 , 35 . D /A 3 

D + L> 3 f — + — - 

18 18 36 144 96 V12 32 A 



ln(A + D) 



A 7 D 7 7 - 

h -D J 

2 2 4 



35 n » 31 „ 
— D + — D 
16 32 



— + -^4- in ( A + ^) 

2D 32A v ; 



44 



487T 2 

A 3 

487T 2 

'A 3 ' 



1 /*q „q\ 3 ~ 21 - 7 „ D /A 

( A 9 - L> 9 ) + D D + D h 

V / inn inn i m 



150 



15 



5 20 8 16 



480 
1 



+ 



V96 80A 
1 



ln(A + D) 



12D 16A 



ln(A + D) 



Qi = 
+ 

e 2 = 



e, 



+ 



A 3 



96tt 2 
"A*" 
D 

1280 ~~ 1^192 
96tt 2 r 1 
~ .90 



— (A 11 - D 
900 V 



ii 



11 
1800' 



-D' 



11 - 77 - L> 3 

Z) + D 

800 4800 640 



A 

— + 



5 



768 A / 
11 



ln(A + D) 
11 



— (D 11 - A 11 ) - —-D 9 + —D 
90 V 180 80 

ln(A + D) 



(A 35 
D+1 — + 



A 5 

768^ + V32 ' 768A 
96tt 2 rl 

~ .6 

7 1 



77 

480' 



_37 

384' 



D 3 — 



A 5 
151 

256 



D + 



(A 11 - D 11 ) + —D 9 - -D 7 + —D 5 - —D 3 + 
V / 12 16 32 128 

315 



24£> 24£> 3 768A 



ln(A + D) 



28 



—1 



—2 — 



327T 2 

A 5 
2 

35D 

327T 2 

A 5 
2 



29 
180 



(A 9 - D ! 



13 3 
315£> 3 + 64 



29 n7 203 n5 20497 . 12 n 

H D 7 D 5 H D 3 D— 

40 160 20160 35 



A + ^ ln(A + D) 



90 V / 9n sn 



1921 



+ 



13 



—3 



35£> ' 315L> 3 
i2n 2 T233 
.180 



+ 



80 
3 



D J 



'A 3 3 . 

1 A + 

24 32 64A 



1680 
ln(A + D) 



2339 
6720" 



D- 



+ 



—4 — 



13 



A 5 
2 

5D ~ 45Z) 3 
32vr 2 T47 



D 



_ 233 7 
40 



1631^ _ 23767^3 



15 . 21 

—A H - 

64 64A 



160 
ln(A + D) 



2880 



449 „ 

D+ 

160 



+ 



—5 — 



A 5 
13 



45 



D 9 - A 
13 21 



47 - 329 . 

— D + D 

10 40 



9611 ~ 663^ 
-D + D+ 



1440 



320 



ln(A + D) 



20D 45L> 3 64A 

A 5 * ^ ' I () 

117 



+ 



10 

+ ttt^t In (A + D) 



3591 n5 11691 n3 



40 



160 



7407 
l320" 



D- 



20D ' 20L> 3 ' 64A 
The tensors which are not integrated in the elementary functions read 



dy 
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where F = y 2 + 1. The coefficients of the obtained polynomials in A and D should be rather pecu- 
liarly organized. Because at large distances the corresponding integrals are the decreasing functions. 
Therefore the coefficients at high powers should be tuned in such a way to cancel each other, up to 
the order necessary for functions decreasing, at expanding D = \J A 2 + 1, when A — > oo. At small 
A the singularity of terms related to ~ ^ what means the sum of all the coefficients, excepting 
the terms proportional to ~ A 3 In A, A In A and the powers of A, should equal to zero. 



References 

[1] C.G. Callan, R. Dashen, and D.J. Gross, Phys. Rev. D 17, 2717 (1978); 
C.G. Callan, R. Dashen, and D.J. Gross, Phys. Lett. B 66, 375 (1977). 

[2] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B 422, 247 (1998); 

R. Rapp, T. Shafer, E.V. Shuryak and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998). 

[3] R. Rapp, T. Schafer, E.V. Shuryak, and M. Velkovsky, Ann. Phys. (N.Y.) 280, 35 (2000); 
G.W. Carter and D.I. Diakonov, Phys. Rew. D 60, 016004 (1999). 



29 



[4] S.V. Molodtsov, A.M. Snigirev, and G.M. Zinovjev, Nucl. Phys. B 106& 107, 281 (2002); 
S.V. Molodtsov, G.M. Zinovjev, Nucl. Phys. B 106& 107, 471 (2002); 
G.M. Zinovjev, S.V. Molodtsov, Yad. Fiz. 66, 1000 (2003), ibid 66, 1389 (2003). 

[5] D. FSrster, Phys. Lett. B 66, 279 (1977); 

E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. B 364, 255 (1991). 

[6] LI. Balitsky and A.V. Yung, Phys. Lett. B 168, 113 (1986). 

[7] S.V. Molodtsov, A.M. Snigirev, and G.M. Zinovjev, Phys. Rev. D 60, 056006 (1999); 

S.V. Molodtsov, A.M. Snigirev, and G.M. Zinovjev, in Lattice Fermions and Structure of the 
Vacuum, edited by V. Mitrjushkin and G. Schierholz, (Kluwer Academic Publisher, Dordrecht, 
The Netherlands, 2000) p. 307; 

G.M. Zinovjev, S.V. Molodtsov, A.M. Snigirev, Yad. Fiz. 63, 975 (2000). 
[8] P. Sikivie, N. Weiss, Phys. Rev. D 20, 487 (1979). 

[9] V.V. Goloviznin, S.V. Molodtsov, A.M. Snigirev, Yad. Fiz. 56, 123 (1993). 
[10] J. Ishida and A. Hosoya, Prog. Theor. Phys. 62, 544 (1979). 

[11] A.M. Polyakov, Gauge Fields and Strings (Contemporary Concepts in Physics: v. 3), Harwood 
Academic Publishers, Chur, Switzerland, 1987. 

[12] D.I. Diakonov and V.Yu. Petrov, Nucl. Phys. B 245, 259 (1984); 

D.I. Diakonov and V.Yu. Petrov, in Hadronic Matter under Extreme Conditions, 
ed. by V. Shelest and G. Zinovjev (Naukova Dumka, Kiev, 1986), p. 192; 

D. I. Diakonov, V.Yu. Petrov, and P. V. Pobylitsa, Phys. Lett. B 226, 471 (1989); 

E. -M. Ilgenfritz and M. Muller-Preussker, Nucl. Phys. B 184, 443 (1981); 
T. Schafer and E.V. Shuryak, Rev. Mod. Phys. 70, 323 (1998). 

[13] D.I. Diakonov and V.Yu. Petrov, ZhETP, 86, 25 (1984). 

[14] H. Aoyama and H. Kikuchi, Nucl. Phys. B 369, 219 (1992). 

[15] J. Smit and A. van der Sijs, Nucl. Phys. B 422, 349 (1994); 

E.-M. Ilgenfritz, S.V. Molodtsov, M. Muller-Preusker, and A.I. Veselov, 
Nucl. Phys. (Proc. Suppl.) B 73, 581 (1999). 

[16] M.I. Polikarpov and A.I. Veselov, Nucl. Phys. B 297, 34 (1988), 
M.-C. Chu, J.M. Grandy, S. Huang, and J.W. Negele, 
Phys. Rev. D 49, 6039 (1994); 

C. Michael, PS. Spencer, Phys. Rev. D 52, 4691 (1995); 
M.G. Perez, T.G. Kovacs, and P. van Baal, 
Phys. Lett. B 472, 295 (2000). 

[17] J.W. Negele, Nucl. Phys. Proc. Suppl. 73, 92 (1999). 



30 



